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Abstract 



We investigate isometric immersions of disks with constant negative curva- 
ture into R 3 , and the minimizers for the bending energy, i.e. the L 2 norm 
of the principal curvatures over the class of W 2 ' 2 isometric immersions. We 
show the existence of smooth immersions of arbitrarily large geodesic balls 

^ ■ in M 2 into R 3 . In elucidating the connection between these immersions and 

the non-existence/singularity results of Hilbert and Amsler, we obtain a lower 
bound for the L°° norm of the principal curvatures for such smooth isometric 
immersions. We also construct piecewise smooth isometric immersions that 
have a periodic profile, are globally W 2,2 , and numerically have lower bend- 

*C\ ■ ing energy than their smooth counterparts. The number of periods in these 

configurations is set by the condition that the principal curvatures of the 
surface remain finite and grow approximately exponentially with the radius 
of the disc. We discuss the implications of our results on recent experiments 



^1" 

Q ■ on the mechanics of non-Euclidean plates. 



Keywords: nonlinear elasticity of thin objects, geometry of hyperbolic 
surfaces, pattern formation, morphogenesis in soft tissue. 



1. Introduction 



The differential growth of thin elastic sheets can generate highly non- 
trivial configurations such as the multiple generation of waves along the edge 
of a torn elastic sheet [l| and the wrinkled patterns on leaves [2|. Recently, 
environmentally responsive gel disks of uniform thickness have been created 
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in the laboratory that mimic this type of growth by differentially shrinking 
in the radial direction when thermally activated in a hot water bath [3| . The 
rate of shrinkage can be controlled, allowing for the creation of regions in the 
disk that shrink at different rates upon activation. Rippling occurs in regions 
that shrink slower than the center of the disk while regions that shrink faster 
obtain a spherical dome like shape. Geometrically, these regions have positive 
and negative Gaussian curvature respectively and by tuning the composition 
of the gel, different curvatures can be selected. This type of growth can be 
used to generate a multitude of different shapes of varying geometries, but a 
mathematical framework that can accurately predict the equilibrium shapes 
is lacking. 

Classically, thin elastic plates have been modeled by the Foppl - von 
Karman (FvK) equations that reduce the full three dimensional equations of 
elasticity to a coupled system of partial differential equations defined on the 
mid-surface of the plate J4J, |5j. These equations are asymptotically valid in 
the small thickness regime if the thickness scales in a particular manner with 
external loadings [6j, and different scalings lead to a hierarchy of different 
plate models that includes FvK |7|. The FvK equations can be obtained as 
the variation of an energy functional. This energy is the sum of stretching 
and bending terms such that the ratio of bending to stretching energy scales 
like the thickness squared |4|. Consequently, thin plates energetically prefer 
to bend out of the plane rather than stretch - a phenomenon easily observed 
in stretched sheets J8( - and this could be the mechanism used by growing 
sheets to buckle out of the plane. 

Non-uniform growth is difficult to model using the FvK equations since 
it commonly generates structures with residual stresses [9| and thus it is 
unclear what is an appropriate reference configuration from which strains 
can be measured. One technique for defining a reference configuration is 
to specify a two dimensional "target metric" gon the mid-surface of the 



sheet from which strains are measured [10|, [111, [12|, |13j. The metric g is a 
geometric quantity that locally determines all of the intrinsic properties of 
an abstract surface such as its Gaussian curvature and geodesies [l4(. Thus, 
we are naturally led to consider the problem of finding isometric immersions 
of g into M 3 [15]. But, the problem is that in general isometric immersions of 
g may not exist or be unique and this model does not incorporate important 
three dimensional effects such as a mechanism for penalizing highly bent 
configurations. 

Instead of modeling non-uniform growth as isometric immersions, a "non- 



Euclidean" elastic energy functional has been recently introduced that incor- 
porates g and reduces to the FvK functional when g is the identity matrix 
[161 ] . This functional has the same energy scaling as FvK with the stretching 
energy measuring L 2 deviations of the immersion from g. As with the FvK 
functional, if the thickness is sufficiently small the stretching energy domi- 
nates and we expect the elastic sheet to bend out of the plane to match the 
metric. In particular, if W 2,2 isometric immersions exist then in the limit 
of vanishing thickness minimizers of the non-Euclidean FvK functional con- 
verge to a minimizer of the bending energy over the class of W 2,2 isometric 



immersions of g [17| . Thus in this model the particular isometric immersion 
is selected by the bending energy. 

An alternative is to decompose the deformation of the sheet into a growth 
tensor and an elastic response, and define the elastic energy in these terms. 
Lewicka, Mahadevan and Pakzad [18| have shown that the two approaches 
are essentially equivalent in the limit that the thickness of the sheet vanishes. 

In this paper we study isometric immersions of g in the disk geometry 
when g specifies a constant negative Gaussian curvature K . Experimentally, 
disks of constant negative curvature obtain an equilibrium configuration with 
a periodic profile of one wavelength with the number of nodes increasing with 



decreasing thickness [19|. Consequently, we will look for low bending energy 
test functions that match the target metric and have a periodic profile. 

Finding such test functions corresponds to the classical problem of finding 
isometric immersions of the hyperbolic plane H 2 into M 3 . Hilbert [20| proved 
that there are no real analytic isometric immersions of the entire hyperbolic 
plane into M. 3 , but this does not exclude the possibility of the existence of 
isometric immersions with weaker smoothness criteria. In fact, the Nash- 
Kuiper theorem states that through the technique of convex integration there 



exists C 1 isometric immersions of H 2 [2l|, |22|. But, such an immersion would 
be too rough since the bending energy requires the surface to be at least 
weakly second differentiable. Furthermore, by Efimov's theorem, we know 



that there are no C 2 isometric immersions of H 2 23j. The admissible surfaces 
with finite elastic FvK energy are W 2,2 and this motivates the search for W 2,2 
isometric immersions which lie between C l and C 2 in smoothness. 

In general finding an isometric immersion of an arbitrary metric is a non- 



trivial problem. But, by Minding's theorem [241 ]. which states that any two 
abstract surfaces having the same constant Gaussian curvature are locally 
isometric, we can cut out subsets of surfaces of constant negative curvature 
which will be isometric to pieces of H 2 . In particular, if S is a surface of 



constant negative curvature and p G S then U = {q G S : d(p, q) < R} is 
isometric to a disk of radius R in HI 2 , provided £/ is smooth. Furthermore, 
it was proven by Amsler that the singularities where S fails to be smooth 
form curves on S 25|. These curves form a boundary beyond which no disks 



locally isometric to H 2 can be cut out of the surface. 

This paper is organized as follows. In Sec. [2] we present the mathemat- 
ical model for thin elastic sheets with a specified metric g in the form of a 
variational problem with a FvK type functional. We show that for isometric 
immersions the functional reduces to the Willmore functional and thus in 
the vanishing thickness limit if an isometric immersion exists the minimizing 
configurations should converge to a minimizer of the Willmore energy over 
the class of W 2 ' 2 isometric immersions 17 . 



In Sec. [3] we review some of the theory of hyperbolic surfaces that will 
allow us to construct isometric immersions in R 3 . In particular, we make use 
of Chebyshev nets (C-nets) which are local parameterizations of a surface 
whose coordinate curves form quadrilaterals on the surface with opposite 
sides having equal length [14J . Originally, Chebyshev used C-nets to describe 
the arrangement of fibers in clothing and he showed that a sphere must be 
covered by at least two cuts of cloth [26| (see [27| for a modern translation). 
As evidenced by the sphere, C-nets in general are only local parameterizations 
of a surface and a global C-net exists if the integral of the mag nitude of the 



Gaussian curvature over a quadrilateral is less than 2tt [28j|. C-nets are 



useful parameterizations in that they lead to a natural way to discretize a 
surface [29J and projections of C-nets onto paper retain the three dimensional 
properties of the surface 30(. Remarkably, for surfaces with constant negative 
curvature, every sufficiently smooth isometric immersion generates a C-net 
on the surface obtained by taking the asymptotic curves as the coordinate 
curves and conversely every C-net on a surface of constant negative curvature 
generates an immersion whose asymptotic curves are the coordinate lines 



14j . The extrinsic and intrinsic geometry of the particular immersion is 
entirely determined by the angle <fi between the asymptotic curves and a 
parameterization is a C-net if and only if is a solution of the sine-Gordon 
equation [14j. Thus, solutions of the sine-Gordon equation generate surfaces 
of negative curvature and the corresponding immersion will fail to be smooth 
when is an integer multiple of tc. 

In Sec. H]we show that isometric immersions of geodesic discs in the hyper- 
bolic plane can by constructed by cutting out subsets of a classical hyperbolic 



surface, the pseudosphere. This is contrary to the claims in |3l| where the ex 



perimentally observed buckling and refinement in wavelength with decreasing 
thickness is attributed to the non-existence of smooth isometric immersions. 
These configurations have very large bending energy even away from the 
singular curve and do not have the rotational n-fold symmetry observed in 



experiments 191 ] . 



A natural question is the connection between this existence result and the 



non-existence result of Holmgren [32[ stating that given a "local" immersion 
of an open set in H 2 into M 3 , there is a maximal domain of finite diame- 
ter on which this immersion can be extended to a smooth global immersion 
with one of the principal curvatures diverging as we approach the boundary. 
In Sec. El we obtain a quantitative version of this result, by showing nu- 
merically that for isometric immersions of a Hyperbolic disk of radius R the 
(nondimensional) maximum principal curvature of the immersion is bounded 
from below by a bound which grows exponentially in the (nondimensional) 
radius. We also identify the surfaces which realize this bound, and they turn 
out to be disks cut from Hyperboloids of revolution with constant negative 
curvature. 

In Sec. El we use the small slopes approximation to find approximate 
isometric immersions with a periodic profile. The small-slopes approximation 



has been widely used to model torn sheets and growing tissues [13|, |33|, |34 
In this setting the problem of finding immersions of g reduces to solving 
a hyperbolic Monge- Ampere equation. The simplest solutions are a family 
of quadratic surfaces with two asymptotic lines intersecting at the origin. 
We then periodically extend the surface bounded between the two lines to 
construct periodic shapes that have the same scaling of amplitude versus 



period observed experimentally [191 ] . The two wave configuration is the global 
minimum of the bending energy over the class of isometric immersions and 
we prove that with decreasing thickness the configurations should converge 
to this shape. 

In Sec. [7] we extend the surfaces constructed from solutions to the Monge- 
Ampere equation by studying hyperbolic surfaces with two asymptotic lines. 
These surfaces are called Amsler surfaces, named after Henri Amsler who 
studied them in his paper on the nature of singularities of hyperbolic sur- 



faces |25|, |35j. Amsler surfaces are generated by a similarity transformation 



of the sine-Gordon equation which transforms it to a Painleve III equation 



in trigonometric form [25|, |35| . From these surfaces, we can create periodic 
shapes as we did with the small slopes approximations that are energeti- 
cally preferred over corresponding disks cut from the pseudosphere and hy- 



perboloids. The bending energy of these surfaces is concentrated in small 
regions near the singular curves and the lines we are taking the periodic ex- 
tension. Furthermore for each n > 2, there is a radius R n ~ log(n) such that 
the n-periodic Amsler surfaces only exist for a radius < R < R n . This 
gives a natural geometric mechanism for the refinement of the wavelength of 
the buckling pattern with increasing radius of the disk. Finally, we present 
a concluding discussion in Sec. El 

2. Model 

We model the stress free configuration of a thin elastic disk of thick- 
ness t and radius R that has undergone differential growth as an abstract 
Riemannian manifold D with metric g. More specifically, if we choose 
a coordinate system for D so that lies at the center of the disk then 
D = {p : d(p, 0) < R}, where d is the distance function induced by g. 
A configuration of the sheet in three dimensions is then given by a mapping 
x : D — > M. 3 , where U = x(£)), the image of the mapping, is taken to be 
the mid-surface of a sheet. The equilibrium configurations of the sheet are 
extrema for the elastic energy functional: 

S[x]= [ \\Dx T -Dx-g\\ 2 dA + t 2 [ \\AH 2 - K\\ 2 dA, (2.1) 

Jd Jd 

where H is the mean curvature of the surface U, K is the Gaussian curvature, 
Dx is the Jacobian matrix of partial derivatives of x, and dA is the area 
element on D [3]. This functional is an asymptotic expression for a full 



three dimensional elastic energy functional per unit thickness [17J. The first 
integral measures the stretching energy and vanishes when x is an isometric 
immersion while the second integral measures bending energy and vanishes 
when the image of x is a plane. 

This functional is a FvK type since it is the sum of stretching and bending 
energies such that ratio of bending to stretching energy scales like t 2 J6[ . If g is 
not Euclidean, i.e. if the Gaussian curvature corresponding to g is non-zero, 
then both of these terms cannot simultaneously vanish in three dimensions 
and consequently the equilibrium configuration will not be able to reduce all 
of its elastic energy even in the absence of external forces. 

In this paper we will focus on (12. ip when g generates a constant negative 
Gaussian curvature K Q < 0. K is the target Gaussian curvature, and we 
can always nondimensionalize such that K = — 1, i.e., we chose the units for 



length as (— Kq)~ 1 ^ 2 . If the surface is (initially) parameterized by u,v then 
u' = y/—Kou, v' = \J—Kqv, are the dimensionless parameterization variables 
and R' = y/—K R is the dimensionless radius. Therefore, if we can find a 
configuration x such that Dx^Dx — g = then, the Gaussian curvature of 
the immersion K = K = — 1 and (12. ip reduces to the following equivalent 
elastic energy functional 



B[x] = / [k{ + k 2 2 ) dA, (2.2) 

Jd 

where k\ and k 2 are the principal curvatures of U. Since k\k 2 = K = — 1, 



this energy is equivalent to the Willmore energy [36[. We will refer to B 
as the (normalized) bending energy of the configuration x, although, strictly 
speaking, the bending energy is the quantity t 2 B and not the quantity B 
itself. 

Since we are integrating over a geodesic circle it is convenient to work in 
geodesic polar coordinates, which are the natural analogs of the radial and 



polar angle coordinates on a Riemannian manifold |14| . In these coordinates 
U is parameterized by the coordinates (r, \l/) which map to a point q 6 U 
satisfying d(q,p) = r and q lies on a geodesic whose initial velocity vector 
makes an angle \l/ with respect to a basis of the tangent plane at p. In this 
coordinate system the metric g has components 

Pn = 1, #12 = #21 = 0, and #22 = sinh(r), (2.3) 

and equation (12. 2 j) takes the form: 

B[x] = / / sinh(r) (k 2 + k%) dr dV . (2.4) 

Jo Jo 



Lewicka and Pakzad |17| have shown that, for immersions x in the space 

W 2 > 2 (D,R 3 ), 

9 „ r . ffifxl, Dx T ■ Dx = g , N 

r-lim r 2 ^[x] = < L J (2.5) 

*^° l+oo, otherwise 

where Z)x T ■ Dx = g implies that x is an isometric immersion. There- 
fore, in the limit as the thickness vanishes, equation ( 12. 4 p is the appropriate 
functional with the set of admissible functions A being sufficiently regular 
isometric immersions of the disk of radius R in HI 2 . In particular, we need 
B\x\ < oo, so the principal curvatures are L 2 functions on the disk D, and 
we need to consider W 2 ' 2 immersions x : D — > M. 3 . 



3. Theory of Hyperbolic Surfaces 

In this section we show how isometric immersions of D can be cut out of 
a surface S with constant negative Gaussian curvature. Let x : D — > S be a 
parametrization of S and for p G S consider the set U = {q G S, d(p, q) < R}. 
By Minding's theorem [24J as long as U is smooth it can form an isometric 
immersion of D. To construct U we need to calculate the geodesies on S and 
also to determine £>[x] we need to calculate the principal curvatures on S. 

To calculate these quantities recall the following facts and definitions from 
the theory of hyperbolic surfaces (see pgs. 683-692 in 14|). An asymptotic 



curve is a curve in S whose normal curvature vanishes everywhere. Moreover, 
since the principal curvatures ki and ki of a hyperbolic surface have opposite 
signs it follows that through each point p G U there are exactly two asymp- 
totic curves that intersect at p. Consequently, we can locally parameterize 
any hyperbolic surface by its asymptotic lines. Such a parameterization is a 
C-net and has a metric g' with the following components 

011 = 922 = 1 aIld 921 = 921 = cos(0), (3.1) 

where <fi is the angle between the asymptotic curves at each point in S. 
Therefore, all of the intrinsic properties of the surface are expressible in 
terms of <f), and for this reason <p is called the generating angle. 

The condition that there must be two distinct asymptotic lines at each 
point in U implies that <ft must satisfy the constraint 

< <f>(u, v) < 7T. (3.2) 



Furthermore, it follows from Brioschi's formula |14[ that must satisfy the 
sine-Gordon equation 

H = -KM*), (3-3) 

and thus solutions to the sine-Gordon satisfying (13. 2 \ generate hyperbolic 
surfaces. In fact Hilbert proved that there is no smooth immersion of H 2 by 
showing that there is no smooth solution to (13.31) that satisfies (13.21) in the 
entire u — v plane. The points on the surface where <fi = nn are precisely the 
boundaries where the immersion fails to be smooth. 

To determine the geodesies on S we need to calculate the Christoffel 
symbols. They are intrinsic quantities determined by the metric and hence 



by 0: 

rii = cot(0& r? 1 = -csc(0)g 

r| 2 = 0, Y\ 2 = 0, (3.4) 

r 22 = -csc(0)§£ r 2 2 = cot(0)§£, 

which gives us the following geodesic equations: 



iir 2 



+ (f) cot(0)^-(§) csc(0)^ = O, 
"(f) csc(0)g+(f) 2 cot(^ = O. 



(3.5) 



The extrinsic quantities of 5 are also expressible in terms of 0. The 
components of the second fundamental form h are given by 

h\i = h 22 = and h\ 2 = h 21 = ±sin(0) (3.6) 

and the principal curvatures are 

k\ = tan 2 (0/2) and k\ = cot 2 (0/2). (3.7) 

Consequently, since U C S is an isometric immersion of D we have that 
equation (12. 4p can be expressed entirely in terms of 0: 

B[x]= / sinh(r)(tan 2 (0/2) + cot 2 (0/2))rfrrf^. (3.8) 

Jo Jo 

We conclude this section by pointing out that a closely related param- 
eterization of S, that is sometimes easier to make calculations with, is a 
parameterization by curves of vanishing curvature or a C-net of the second 
kind [14j. If x is a C-net, then a C-net of the second kind y is given by 
y(Vi£) = x {J21 2~ )• ^ ne com P onen ts of the metric g' corresponding to y 
are 

0n = cos 2 (6>), g' 12 = g' 21 = 0, g' 22 = sin 2 (6»). (3.9) 

Since g[ 2 = g' 21 = 0, the parameter curves of this parameterization are or- 
thogonal and in fact they bisect the asymptotic curves. Therefore, 9 = 0/2 
and thus we can determine the generating angle by such a parameterization. 

4. Elastic Energy of the Pseudosphere 

To illustrate the process of calculating the elastic energy of an isometric 
immersion of D we consider a commonly known hyperbolic surface, the pseu- 
dosphere. The pseudosphere is given by the following C-net of the second 



kind |34 

, .. / cos(0 sin(0 . . A .. „ N 

\cosh(77) cosn(77j J 

The mapping is singular on the curve rj = 0, which is where the surface 
fails to be a smooth immersion, and consequently we will only consider the 
"upper half 77 > of the pseudusphere. Furthermore, since the sheets we 
are modeling have the topology of a plane and not a tubular surface we will 
think of the pseudosphere as a multi-sheeted surface. 

Directly calculating, the metric coefficients are given by the following 
simple expressions: 



g' u = tanh 2 (r/), g 22 = ^— and g[ 2 = g' 21 = 0. (4.2) 



1 
cosh 2 (77) 

Consequently, by equation (13.91) the generating angle has the form 

0(77,f) = 4arctan(e" r '), (4.3) 

and thus by ( 13.71) the principal curvatures are 

k\ = tan 2 (2 arctan(e _r? )) = 2— an d k 2 = cot 2 (2arctan(e _r? ) = sinh 2 (?7). 

sinh (rj) 

(4.4) 
Therefore, by equation (I3.8P the elastic energy of a disk U centered at (770, 0) 
and lying in the pseudosphere is simply 

#[yl = / / sinh(r) ( sinh 2 (?7(r, ^)) H *— — )drd^. (4.5) 

Jo Jo V UV JJ sinh 2 (r7(r,^)); l ; 

Now, to determine how 77 depends on r and \l/ and the shape of a disk 
U we need to calculate the arclength of geodesies on the pseudosphere. The 
arclength of a curve £(77) starting at (770,0) and terminating at (77/, £/) is 
given by the functional 



L[£(t])] = / f Wtanh 2 (77) + sinh 2 (r/) ( -£ J dr; (4.6) 

and the geodesies that extremize this functional are given implicitly by 

cosh 2 (r/) + (f + C) 2 = D, (4.7) 

10 



where C and D are determined by the condition that the geodesic passes 
through (770,0) and (»7/,£/). 

Define the function £(77, £) to be arclength of the geodesic starting at 
(770,0) and terminating at (77, £). Since the geodesies fail to be a function 

of 77 at the critical value I 77* = arccosh(vZ>),£* = Wcosh 2 (77 ) — cosh 2 (77j) J 

the arclength of a geodesic will have to be computed over both branches of a 
square root. Assuming that 770 < 77, if 770 > 77 then £(77, £) can be computed 
by switching 770 and 77, we have by equation (14.61) that 



c(v,t) 



P j ^m dt 



VO V D~cosh (t) 



JVO VV D-cosh 2 ^) + V D-cosh 2 (t)J m 



hi 



Mjj) \ 



if 


< e < (r 


if 


e > (d 2 


if 


^ = 



cosh(r;o)y 

(4.8) 
By making the substitution z = sech(t), and taking into account the case 
770 > 77, the integrals can be explicitly evaluated. The contour £(77, £) = R 



forms the boundary of the disk U, (see figure 1(a) ). 



Furthermore, the geodesies on the pseudosphere itself are given by x(t] = 
arccosh (D — (£ + C) 2 ) ,^). The tangent vector jP evaluated at 77 = 
gives us the tangent vector of this curve which we can use to calculate 

the angle a geodesic makes with respect to the basis ( 2 , q ) of the 
tangent plane at the center of the disk by taking inner products. By making 
the substitution D = cosh(77 ) 2 + C 2 , this calculation gives us that the angle 
\l/ satisfies 

cos W = v ^ csch ^°) U 9) 

C ° Sl j l + 2C 2 + cos(2£o)' [ ' 

Consequently, by specifying a value for \1/ we can calculate C and therefore 
use equation (14. 8p to determine 77 for a particular value of r and \1/. This 



allows us to color the disk by bending energy density (see figure 1 (b) ) and 



numerically integrate (14.51) over the disk to determine how the bending energy 



scales with R for various values of m (see figure 1(c) ). 

Now, a disk centered at (/70, 0) cannot be made arbitrarily large since it 
will eventually meet the singular curve 77 = 0. We can prove that the largest 
disk on the pseudosphere with this center has a radius of R m = ln(cosh(77 )). 
To prove this note that the curve £ = is a geodesic, it trivially solves the 
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(a) (b) 

Scaled-Beiiding-Energyvs.'GeodesicR.adms'of'Pseudosphere 




Figure 1: (a) The pseudosphere colored by the geodesic distance from an arbitrary point 
p. Geodesic disks of constant negative curvature can be formed from the contour data 
and correspond to isometric immersions of geodesic disks in the hyperbolic plane. |(b)| 
A representation of the geodesic disks colored by k\ + fcf in which the radius and polar 
angle correspond to the geodesic radius r and polar angle ^. These disks are centered 
at 770 = 1.94,2.54,2.92,3.21 respectively. These values where chosen for comparison with 



the periodic Amsler surfaces in figure ® (c) Plot of the scaled bending energy of these 
disks as a function of the geodesic radius. We can see from figure [(b)] that for moderately 
sized disks the bending energy is large near the top of the disk leading to the very large 
energies. 
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Euler-Lagrange equation, and therefore 



I/O 



£(0,0)=/ tanh(u)dr] = \n(cosh(r] )). (4.10) 

Jo 

Now, let £(77) be another geodesic that terminates at a point (0,77*) on the 
boundary. Then, 




2/ \ , ■ , 2, , / <% 



2 



£(!)/;)= / \/tanh ^(r 7 ) + sinh fa)!/) ^ 

> / tanh(?7) ci?7 = ln(cosh(?7o)). 
Jo 

Consequently, the contour £(77, £) = ln(cosh(r7o)) meets the curve 77 = at 
only one point and thus ln(cosh(77o)) is the radius of the largest possible disk 
centered at (770,0). Therefore, we can see that for arbitrary R we can find 
a value 770 such that D is isometric to U, and thus by varying the center of 
the disk we can create arbitrarily large stretching free configurations. This 
proves the following proposition 

Proposition 4.1. Let D be a disk of radius R in the hyperbolic plane. There 
exists a smooth isometric immersion x : D — > U C M, 3 such that U is a subset 
of the pseudosphere. 

5. Lower Bounds for the Curvature of Smooth Isometric immer- 
sions of Hyperbolic disks 

In this section we explore numerically how the principal curvatures of 
smooth isometric immersions scale with the size of the disk. In particular, we 
show that the maximum principal curvature on the disk grows exponentially 
in ^/—KqR, where K is the target Gaussian curvature. 

Let p denote the center of the disk, and let us choose the asymptotic 
coordinates u and v such that p is the point (0,0). If |w | + l^ol < R, it 
follows easily that there exists p' in D with coordinates (uq, i>o) since d(p, p') < 
\ u o\ + \vo\ < R. Let q, r, s and t denote the vertices of a "square" denoted 
[qrst] in asymptotic coordinates given by the intersections of the asymptotic 
lines u = -R/(2 + v),v = -R/(2 + v),u = R/{2 + v),v = R/(2 + u) for 
some given v > 0. 
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We redefine the variables 

, 1 (2 + rj)u , 1 {2 + rj)v A\K \R 2 , 

u — — , if = — , A = - — , d> [u , v ) — 4>{u, v ), 

2 2i? ' 2 2# ' (2 + r/) 2 ' ^ v ' ; vv ' ; ' 

and then drop the primes to transform equation (13.31) into 

4>uv = Asin(0) on the domain (0, 1) x (0, 1) (5.1) 

A smooth isometric immersion corresponds to a smooth function : (0, l) 2 — > 
(0, 7r) satisfying the Sine-Gordon equation (15. ip . 

We want to consider the variational problem for minimizing the LP norms 
of the principal curvatures. The natural mode of converge is then LP con- 
vergence and this does not preserve the class of smooth functions. We thus 
enlarge the admissible set. Let Ai denote the set of all measurable functions 
on the unit square with values in [0, ir] let A C Ai be the set of measur- 
able functions <p : [0, l] 2 — > [0,7r] that are distributional solutions of the 
Sine-Gordon equation, i. e 

cj)g uv dudv = A / sm((J))gdudv, \/g G C^°((0, l) 2 ). 
[0,1] 2 J[o,i] 2 

Clearly A C U>([0, l] 2 ) for all 1 < p < oo. 

Remark 5.1. If <p n is a sequence in A and n — >■ in LP, the mean value 
theorem implies that || sin(0 n ) — sin(0)|| p < ||0„ — 0|| p —¥ ; so it follows that 
4> G A. Consequently, A is a closed subset of L p ([0, l] 2 ) for all 1 < p < oo. 

Remark 5.2. (Compactness) Since M. consists of bounded measurable func- 
tions, for every sequence 4> n G M., there exists a subsequence (p nk and a limit 
point <fi* G Ai such that 4> nk — ^ <p* weakly in all L p ([0, l]) 2 for 1 < p < oo 
and <p nk — *• <fi* (converges weak-*) in L°°([0, l] 2 ). Note that the limit <p* is 
independent of p. In fact, on M., the weak L p topologies for 1 < p < oo and 
the weak-* L°° topoloy are all the same \3nj . 

We are interested in a lower bound for the principal curvatures of a Cheby- 
shev patch, i.e., a portion of a surface of constant negative curvature covered 
by a single coordinate system given by a Chebyshev net. We thus define 

M(A) = inf sup max 

<t>^ A a;e[0,l]x[0,l] 
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cot I — I , tan — 
1 2/ ' \2 



While we can estimate M(X) by direct numerical computation, we use a 
slightly different approach that yields some more information. Since cot (0/2) + 
tan(0/2) = 2csc(0), we get k\ + k\ = 4csc 2 (0) - 2 = 4cot 2 (0) + 2. This 
motivates the definition 

J p (A) = inf ||cot 2 (0)|| p 

where ||.|| p denotes the LP norm. This corresponds to the variational problem 
for functionals F p : M. — > [0, oo] given by 

; , (n] = i M cot 2 0||p, (peA. 

+oo, otherwise. 

Note that M(A) and I P (X) are for a unit square in asymptotic coordinates, 
and this is not the same thing as a geodesic disk, although the two can be 
related. 

For all p < oo, the map / i— > \\f\\ p is differentiable at / ^ 0. Also, for any 
function / : [0, 1] x [0, 1] — ¥ R, the p-norm ||/|| p is a nondecreasing function 
of p. It is easy to see that, as a consequence, I p is non- decreasing in p. 

A natural question is whether there exists a minimizer for F p . In this 
context we have 

Proposition 5.3. Let <p n be a minimizing sequence for F p , i.e., <p n G A and 

F p (4> n ) -^ I P - If 4>n (or a subsequence) converges pointwise to <fi* , <fi* is a 
minimizer for F p . 

Proof. If (j) n (or a further subsequence) converges pointwise, it follows from 
dominated convergence that \\<f> n — <f>*\\ p — > 0. Since A is closed in L p , it 
follows that 0* G A. Also, by Fatou's lemma, we have 

L p <F p {<j ) *)<\imintF p {<P n ) = L p 

showing that the minimum is attained. □ 

Remark 5.4. The function 9 \- > — Asin(#) is convex on (0, 7r). Therefore, 
the map 

M 3 i-> -A / sm{<f))hdudv, h G C c °°((0, l) 2 ), h > 0, 
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is weakly lower semi- continuous. If <f> n G A is any sequence and n — *■ 0* 
weakly in some L p , 1 < p < oo or converges weak-* in L°° , <p n —*■ 0* it follows 
that 

/ <p*h uv dudv - A / sin((j)*)hdudv < 0, V/i G C c °°((0, l) 2 ), h>0. 

J[0,l] 2 J[0,1] 2 

So, in general, 0* ^ A, and thus, the functional F p are not lower semi- 
continuous with respect to weak or weak-* convergence. Thus, it is not enough 
to have weak convergence of the minimizing sequence for the existence of a 
minimizer. The pointwise convergence in the proposition is also necessary 
for using the direct method to show the existence of a minimizer for F p . 

The previous remark suggests considering the problem of minimizing the 
relaxed energy F p , namely the largest weakly lower semi-continuous function 
less than or equal to F p on M. 38|, |39[ . Let B denote the weak closure of A 
in some LP (the precise value of p' is irrelevant as indicated in Remark I5.2J) . 
The relaxation of F p is defined by 

~ _ \ inf{liminf n F p (0 n ) : n ->• 0} e B. 

l+oo, otherwise 

where the infimum is over all sequences which converge weakly to 0. Clearly, 
A C B C M. and F p is weakly lower semi-continuous by construction. It 
thus follows from Remark 15.21 that we can prove the existence of a minimizer 
for F p by the direct method in the calculus of variations. 

The question of the existence of a minimizer for F p can thus be posed 
in terms of the nature of the minimizer of F p which always exists. If there 
exists 0* G A which is a minimizer for Fp, then 0* is also a minimizer for F p . 
Conversely, if all the minimizers of F p are in B \ A, it follows that F p does not 
have a minimizer. Studying this question involves computing the relaxation 
F p which in turn needs a careful analysis of the potential oscillations in 
weakly convergent sequences of solutions of the Sine-Gordon equation. We 
are investigating these issues in ongoing work. 

5.1. A numerical investigation of lower bounds for the curvature 

We can discretize the Sine-Gordon equation on a iV x N grid covering 
the unit square, and minimize the V norms of cot 2 over all solutions of this 
discrete Sine-Gordon equation. This is a finite dimensional problem with a 
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coercive energy, so there always exists a minimizer on the N x N grid. 



Fig. 2(a) displays the numerically obtained results for this minimization. 
Observe the lack of high frequency oscillations in the numerically obtained 
minimizers. In fact, increasing the mesh size suggests that as N — > oo, the 
numerically obtained minimizers converge pointwise on [0, l] 2 . If this were 
indeed the case, Proposition 15.31 implies the existence of minimizers for the 
functionals F p on the admissible set A. 

We will henceforth assume that the infimum is attained, i.e., for any 
given p and A, there is an admissible function P (A) = P G A such that 
I p = || cot 2 (0 p )||p, it follows that 

Ip<Ioo< || cot 2 (0p)||oo. 

Thus, numerically determining the function <fi p which minimizes the p-norm 
of cot 2 (0) in the admissible set, gives both lower and upper bounds for i^, 
and the difference between these bounds gives an estimate for the error in a 
numerical determination of 1^. 

Once we determine loo (A), we can compute M(A) as the following argu- 
ment shows: 



max(|fci|, \k 2 \) 



\h\ + I&2I \\h\ - \h 



k\ + 2\k x k 2 \ + k\ k\ - 2\k x k 2 \ + k 2 2 

4 + V 4 



rCi ~\~ Zi ~\~ /Co / rui Zi ~]~ /Co 

4 4 



v/cot 2 (0) + 1 + v /cot 2 (0) = s(cot 2 (0)) 



where g(x) = \fx+ yl + x is a monotone function with a monotone inverse. 
Consequently, 

M(A) = g(lM)) = v / ^UAj+ v/^A) + 1, 

and the value of M(A) is attained on (p^, the extremizer for 1^. 

The same argument can be applied to any function ((<ft) with the property 
that there is a continuous monotone function h with a continuous monotone 
inverse such that h(((<p)) = cot 2 (</>). In particular, if C(0) = {<p — 7r /2) 2 , there 
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is such a function h, and it follows that 



are min max cot 2 (d>) 
<t>&A x&D 



argminmaxmax[tan(0/2), cot(0/2)] 
4>^A xeD 



arg mm max 



TT 



i2 



(5.2) 



We can solve any of these minimax problems to determine the extremizing 
function (p,^. For our numerics, we use the third formulation of the problem, 
i.e., finding the minimax for the square of the deviation of <p from tt/2. Also, 
we note that, 1^ directly estimates the maximum bending energy density 
k\ + &2 and I\ estimates the total bending energy. 



Figure 2(a) shows the numerically computed values of I p and || cot (<p p 



as a function of p = 2 n , nGN, for A = 4. Figure 2(b) and figure 2(c) shows 



the numerically determined extremizers <p p as well as the quantity cot 2 (0 p ) 
for A = 4 and p — 16. The figures for other values of A are similar in char- 
acter. The key observation from these figures is that as p gets large (p p (u, v) 
numerically approaches a limit that is purely a function of the combination 



u + v (see figure 2(d) ). This motivates the ansatz (poo(u, v) = ip(u + v). 

Substituting this ansatz in equation (15.1 ft yields ip" = Asin(-?/>), which is 
the Hamiltonian motion of a unit mass in a potential V(ip) = \cos(ip). The 
trajectories in phase space (ip, ip') are given by conservation of the energy, 
ip' 2 /2 + Acos(^) = E (see figure I5.ip . Note that there are three distinct 



types of surfaces of revolution of constant negative curvature [40[ and for 
each of these surfaces <p is a function of u + v only. Therefore, the three dis- 
tinct trajectories in the phase plane correspond to the three types of surfaces 
of revolution. The pseudosphere corresponds to the separatrix, hyperboloid 
surfaces correspond to closed orbits, and conical surfaces correspond to un- 
bounded orbits. 

We seek solutions on the unit square, 0<t = u + v<2, that are 
bounded away from ip = n on the interval [0,2]. If we want trajectories 
that minimize the maximum value of cot 2 (^), it is clear from the symmetry 
of cot(ip) we are interested in the closed orbits. Furthermore, we seek a 
piece of these trajectories that satisfies the further symmetry requirement 
ip(0) = ip(2) = Ti — 6,-0(1) = e, and ^'(1) = since if these conditions are 
not met we could select a different piece of the trajectory that could further 
minimize cot 2 (0). This trajectory is depicted in bold in figure loTTl 



Estimates of Too (4) 




(a) 





(b) 



(c) 




(d) 



Figure 2: (a) Plot of I<in and || cot (</>2™)|| as a function of n with A = 4. The collapse of 



these data points onto a single line confirms that these two quantities are valid approxima- 
tions to loo for large enough n. |(b)| Plot of 4>ie as a function of u and v with the contour 
plot of 016 underneath the surface. We can see from the contour plot that it appears that 
lines u + v = constant form the contours suggesting that (p^ is a function of u + v. (c) 



cot (</>i6) plotted as a function of u and v again with the contour plot cot ((pie) plotted 
underneath the surface. |(d)| Plot of p (u, v) as a function of u + v for p = 2, 4, 8, 16 the 
collapse into a single curve again indicates that 0oo is a function of u + v only. 
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Figure 3: Phase portrait of the differential equation ip" — Asin(V'). The three different 
types of orbits correspond to the three different types of surfaces of revolution of constant 
negative curvature. The piece of the orbit outlined in bold is the trajectory ip : [0, 2] — > 
(0,7r) that minimizes the value of cot 2 (V>) over all trajectories in the phase plane defined 
over the same time domain. 
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The value E of a minimizing trajectory can be computed as E = [ip'(l)} 2 /2+ 
\cos(ip(l)) = Acos(e), and the relation between e and A is determined by 
requiring that the time it takes to get from ip = n — e to ip = e is 1, i.e., 

1 (5.3) 



V / 2A(cos(e) - cos(0)) 

We can solve this equation numerically, and thereby determine oo ( , u,t') = 
ip(u + v), provided that our ansatz is valid. 

Generating an asymptotic series for ( |5.3p . we have for small e that 

-log(e)+log(8) + Q(e) _ 

Vx 

Rearranging, we obtain 

cot(e) ~ e- 1 + 0(e) = \ exp[V\] + 0[e-^} (5.4) 

8 

This implies ioo(A) ~ ^exp[2vA] + 0(1). Finally, if we set v = we have 
in terms of \/—KqR that 

I^^KR) ~ ^- A e 2 ^ R + 0(1). (5.5) 

5.2. Elastic energy of hyperboloids of revolution 

We conclude this section by comparing the bending energy of disks cut 
from hyperboloids of revolution with the pseudosphere. The family of hy- 
perboloids are parameterized by the following family of C-nets of the second 
kind 

y(?7, £) = - (dn(?7, b 2 ) cosibv), dn(r], b 2 ) sin(fe^), u — E(am.(r], b 2 ) \b 2 )) , (5.6) 
" b 



where dn, am, sn denote the usual Jacobi elliptic functions [41(, E is the 
elliptic function of the second kind 41] , and < b < 1 14j. The metric and 
principal curvatures are given by 

g' n = b%n 2 ( V , b 2 ), g[ 2 = g' 21 = 0, g' 22 = 1 - &Wfa, b 2 ) (5.7) 

k 2 = 922 l-& 4 sn 2 (77,6 2 ) 2 = b 4 sn 2 (r),b 2 ) 
1 g' u 6 4 sn 2 (r/,6 2 ) ' 2 1 - & 4 sn 2 (??, b 2 )' { ' 
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The Christoffel symbols for this parameterization are 

-pl _ cn(r;,b 2 )dn(??,fc 2 ) p2 _ n 

1 11 _ snln.b 2 ) ' L 11 _ U 



r 1 

1 12 

r 1 

1 22 



cn(r),b 2 )dn( V ,b 2 ) 
sn(r),b 2 ) ' 

o, 

cn(r],b 2 )dn(ri,b 2 ) 
sn(r],b 2 ) 



r 2 

1 12 



-b 4 sn(t ? ,b 2 )cn(r ? ,b 2 )dn(r ? ,fe 2 ) 
l-6 4 sn(r ? ,fe 2 ) ■ 



(5.9) 



ri 2 = o, 



Therefore, the geodesic equations are 

d 2 rj cn(r/, b 2 )dn(i], b 2 ) 
dt 2 sn(?7, b 2 ) 



dtj \dt 



d 2 £ 2b 4 sn(r], 6 2 )cn(7/, 6 2 )dn(r/, 6 2 ) 
= ~dt 2+ l-6 4 sn 2 (r7,6 2 ) 



(5.10) 
(5.11) 



It is easy to see that the principal curvatures diverge when it = am(^, b 2 ). 
Consequently, the center of the disk is located at the coordinate i] Q = F(tt/2, b 2 ), 
F is the elliptic function of the first kind, and thus the relationship between 
e and b is given by 



- = arctan 
2 



l-6 4 sin 2 (am(F(7r/2,6 2 )) 



6 4 sin 2 (am(F(7r/2,6 2 )) 
Solving this equation yields 

b= v/cos(e/2). 



arctan 



1-6 4 



b [ 



(5.12) 



(5.13) 



Therefore, for a particular radius R we can numerically solve equation (15.31) 
to obtain b. Once, we determine b we can then numerically integrate the 
bending energy 



B\y] 



2n pR 



JO 



1 — b 4 sn 2 (r], b 2 ) b 4 sn 2 (r], b 2 ) 



6 4 sn 2 (ry,6 2 ) l-6 4 sn 2 (??,6 2 ) 



sinh(r) drd%>, (5.14) 



following a similar procedure as we did with the pseudosphere. These results 
are summarized in figure 15.21 

6. Small Slopes Approximation 

The configurations formed from the pseudosphere and hyperboloids can 
have very large bending energy concentrated at the center and top edge of 
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-2 -10 12 




-2 -1 



(a) 



Sca]ed'Bendiiig-EnErg>'-VB'Radius'foT'Disks-Cut'FToiT]'Hyperboloids 




(b) 



Figure 4: (a) A representation of the geodesic disks cut out of hyperboloid surfaces 



colored by kf + k%- The radius and polar angle of these disks correspond to the 
geodesic radius r and polar angle '&. These disks each have a maximum radii of 
R = 1.2654, 1.8505, 2.2342, 2.5199. These values were chosen for comparison with fig- 
ures Q] and |(b) | The scaled bending energy of disks cut from hyperboloids plotted versus 
geodesic radius. 
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the disk and these configurations do not have periodic profiles. In order to 



understand the experimental observations for negatively curved disks [19 
we thus need to identify periodic configurations with a lower elastic energy. 
We begin by considering the small slopes approximation. In this approx- 
imation, valid for e = y/—K <C 1, the configuration is assumed to have a 
parameterization x : D — > M. 3 of the form: 

x(w, v) = (u + e 2 f(u, v), v + e 2 g(u, v), eu(u, v)), (6.1) 

where f , q G W 1)2 (D), u G W 2 ' 2 (D) and u, v are Cartesian coordinates on D 
\4, UM, |34J- That is, as a mapping, x = i + ecu + e 2 x, where i is the standard 



immersion IR 2 — > M 3 , u is the "out of plane" deformation, i.e., u maps into 
the orthogonal complement of i(M 2 ) and x = i ° (/jfl 1 ) is the "in-plane" 
displacement. The target metric g takes the form 

ds 2 = g = dr 2 + fr 2 + —\ d^ 2 = (dr 2 + rW) + e ^ r ^l ( 6 .2) 

(V XIV 11 \ 
du 2 + 2—dudv dv 2 J . (6.3) 

That is, g = go + e 2 gi is an asymptotic expansion for the target metric where 
go is the Euclidean metric on M 2 . 

If we define the in-plane strain terms by 

= 2 §/ + (^) 2 -^ 

du \ au I 3 

_ 9f_ _i_ 9g_ _|_ du_dw_ , uv (R A) 

dv ~ r du ~ r du dv "•" 3 V J 

dv \ dv J 3 

then using the above asymptotic expansions in the elastic energy ( 12. ip and 
keeping the lowest order terms, we get 

f [x] = e~ 2 £[x] 

= f (\\Dx T + D X + DoJ T -Du- gl \\ 2 + t 2 \\D 2 uj\\ 2 )dA+0(e 2 ,t ) 
Jb(r) 

(6.5) 



/ (tJj + 2 7l 2 2 + 7 2 2 ) dudv + t 2 \\D 2 u\\ 2 dudv + C(e 2 , t°), 

JB{R) " JB{R) 

(6.6) 
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where B(R) is a (2 dimensional Euclidean) disk of radius R and dA is the 
Euclidean area element. The stretching and normalized bending energies 
under this approximation take the forms 

5[x] = / ( 7l 2 ! + 2 7l 2 2 + 72 2 2 ) dudv (6.7) 

Jb(r) 

58 [x] = / ||Z> 2 u;|| 2 ckcfr;, (6.8) 

with the appropriate elastic energy given by 

£ t [x] = S[x]+i 2 <B[x]. (6.9) 

We are using the symbols S and 03 to indicate that these are the lowest order 
contributions to the stretching and bending energies. 

In this approximation, the Mean and Gaussian curvatures to the lowest 
orders in e are given by 

H = e Ti[D 2 u}, K = e 2 Det[D 2 u] 

so that a necessary and sufficient condition for the surface to have constant 
negative curvature K = — e 2 is that the out of plane displacement u should 
satisfy the following hyperbolic Monge- Ampere equation 

d 2 u d 2 u / d 2 u y _ 1 
9m 2 <9t> 2 \ dudv J 

Once a solution to f)6.10p is chosen we can match the metric g' to the target 
metric g by solving for the in plane stretching terms / and g. 

A class of solutions to f)6.10p is the following one parameter family of 
quadratic surfaces 

u a (u, v) = - ( au 2 v 2 ), (6.11) 



parametrized by a G M. + |42|. These solutions are saddle shaped such that 
oo vanishes on a pair of lines passing through the origin of the u — v plane 
intersecting at an angle 9 = ±arctan(a). 

To create periodic structures we first consider the function u* a (6,r) = 
u a (r, (8 — arctan(a)) which is simply a rotation of u a that aligns the asymp- 
totic lines with the u axis and the line 9 = 2arctan(a). Then, we take the 
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odd periodic extension uj a of u*(9, r) with respect to 9 about the two asymp- 
totic lines. Now, since the period of u a is 4arctan(a), to guarantee that the 
shape is continuous at 9 = we must have that a satisfies arctan(a) = 7r/2n, 
where n is the number of waves in the configurations profile (see figure [5]). 




Figure 5: The one parameter family of solutions u> a are saddles with asymptotic lines that 
vanish along the lines through the origin described by the polar angles 9 = ± arctan(a) 
colored in red. Periodic solutions to the Monge- Ampere equation with n € {2, 3, . . .} peri- 
ods can be constructed by taking odd periodic extensions of the solutions u} ta , n (-K/2n)( u , v ) 
about the line 9 — n/2ir. These solutions are C l {D) and are smooth everywhere except 
on the asymptotic lines. 



Since u a is smooth and vanishes along the two lines we are taking the 
extension about, it is a basic result of Fourier analysis that uJ a has continuous 
first partial derivatives but a jump in the second partial derivative with 
respect to 9. Consequently, these lines are singularities in the sense that 
the surface is not smooth across these lines but the bending energy remains 
continuous and is finite. Furthermore, since the metric only involves first 
derivatives of the map x the stretching energy will still vanish on these lines. 
Therefore, these surfaces are valid "isometric" immersions for the small slopes 
approximation i.e., they are W 2,2 isometric immersions. The elastic energy 
of a disk with radius R and n periods is simply: 



33 57 



<£ 



2n 



|£>V 



\ 2 r dr d9 



ttR f tan 



-) + cot 2 ( f ) 
2nJ \2nJ 



(6.12) 



For any solution of (I6.10p . the product of the eigenvalues of the Hessian 



D u is k\ki = — 1. Consequently, 



\D 2 u\ 



Ki-\ \ n^o — /vi i" fCi ^_ ^ 
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with equality only if k 2 = k\ = 1. With n = 2 (i.e. a = 1) in the periodic 
profiles from above, the Hessian of the surface has eigenvalues k\p = il- 
Consequently, the 2-wave saddle u(u, v) = \{u 2 — v 2 ) gives a global minimum 
for the bending energy over all the solutions of (I6.10p . 

Using this as a test function for the elastic energy in (I6.6p . we get an 
upper bound for the energy of the minimizer which scales like t 2 . It is a 
standard argument to show that for all t > there exist minimizers for the 



reduced elastic energy E t |37j|. Further, if the appropriate r-limit of E t exists 
then these minimizing configurations and their derivatives will converge as 
t — > to a limit configuration which satisfies satisfies f)6.10p and minimizes 



the bending energy (16.81) 43] . We provide a more detailed argument proving 
these statements below. 

6.1. V -convergence and convergence of minimizers 

Recall that on a metric space X a sequence of functionals F t : X — > M. 
T-converge as t — > to a functional F : X — > R, written T — lim t ^ F t = Fq, 
if for all x £ X we have that 

1. (Liminf Inequality) for every sequence Xj converging to x 

F [x] ^liminf-Fj-fo]; (6.13) 

3 

2. (Recovery Sequence) there exists a sequence Xj converging to x such 
that 

F [x] >limsupF t [4 (6.14) 

3 

Proposition 6.1. Let Ft be the sequence of functionals onA = W 1 ' 2 (B(R))x 
W l > 2 (B(FL)) x W 2 ' 2 (B(R)) defined by 

F ( [x]=F f [(/, 5 , U )] = ^[x], (6.15) 

where we are naturally identifying a configuration x with an element (f, g, u) £ 
A. Then 

r T T?\(t W T?\(f M J®[ W ] ^711=712 = 722 = 

r-hmF t [(/ ; 5,a;)] = F [(f,g,u)\ = < . (6.16) 

f ^° oo otherwise 
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Proof. Let x t = (f t , g t , u t ) G A such that x f -^ x = (/, g, u). 

1. It is a basic property of weak convergence that W/ G L 2 (B(R)) if 77^ — ^ 77 
in L 2 (B(R)) then \\r]\\ L 2 < liminf t ^ ||^t||_L2 and 3C > such that 

IMU a < C- 



2. It follows immediately from item 1 that 



03 [x] < liminf<B[x t l. 



3. Since B(i2) has finite area it follows that ^ G W 1 'P(S(i2)) for 1 < 
p < 2 and thus by Rellich's compactness theorem we have that ^ 
converges strongly to |^ in the L q norm for 1 < q < 00. Therefore for 
1 < q < 00 we have by the reverse triangle inequality that 



du t 
du 


Li 


du 
du 


Li 


< 



duh 
du 



du 
du 



Li 



and thus liuu.n ||%M| r „ 

l ^ u II au II Li 

limits and inequalities: 



|§^IL„. Therefore we have the following 

I au II Li c 



lim 



dft 



du 



> 



2 

L 2 - 



L 2 



11 



lim 

i->0 



B(R) 



du 



du t 
du 



dA 



lim 



B(R) 



du 



du t 
du 



\du) 



dA 



,. f dft fduV , A 
+ lim / -^i ( _ cL4 

^JBiR) du \duj 

dj_ (du\ 2 rfA 



111 



lim 



9 MdA 



'-0 / B(fl) 9m 3 



B(fl) 



9 /^ ^ 
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IV 



du t 

du 


4 
L 4 


dec 
du 



lim 

t->o 



B(R) 



lim 

i->0 



Scut \ f 2 



du 



dA 



B(R) 



L 4 



fdu\ 2 v 2 , , 



Therefore by items (i-v) it follows that 



lim inf 
t-*o 



B(R) 





B(R) 



2 df_ fdu\ _ f_ 
du \du J 3 



Identical arguments can be used on the other strain terms to show that 

<S[x] < lim inf S [x t ] . 



If <S[x] = then by item 2 we have that 

liminf Fjx t l = lim inf — <S[x*l + lim inf 03 [x+l > Q3fxl, 
t-to t-+o t 2 *->o 

while if S [x] ^ then by item 3 
liminf i^fxj] > lim inf — «S[x t ] + lim inf 03 [x t ] > lim inf — lim inf 5 [x 4 ] = oo. 



t-+o 



t-+o t 1 



t-^-o 



' t->o t 2 *->o 



Therefore, lim inf t ^ F t [x t ] > F [x] proving property 1 of the definition. Fi- 
nally, if y G A then by selecting the sequence y t = y we prove property 
2. □ 



Lemma 6.2. Let at 6e the sequence in 
limt-xjGt exists and < lim^oQi < 2-k. 



defined by at = min^Ft. Then, 



Proof. Fix s, t > such that s < £ and let x s ,x t e ^4 such that min_4F t 
F t [x t ] and min^F s = F s [x s ]. Then, 



F t [x t ] < F t [x s 



1 

f2 



5[x s ] + 03[x s ] < -U[x s ] + 03 [x s ] = F s [x s 



Therefore, the sequence at = min^ F t is monotone increasing as s — > and 
by (16.12p . with n = 2, satisfies < a t < 27ri?. Therefore by the Bolzano- 
Weierstrass theorem at converges. □ 



dA. 
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Definition 6.3. Define Aq C A by 

A = {(f,g,u)eA: /(0) = g(0) = w(0) = and ^- 



(0,0) 



9a; 



(0,0) 



Remark 6.4. TTie above definition fixes the origin of a configuration and 
aligns the normal at the origin with the z-axis ofM. 3 . By a rigid translation 
and rotation any element x G A can naturally be identified with an element 
ofAo. 

Theorem 6.5. If x. t G Aq is a sequence such that min_4 F t = F t [x. t ] then 

2-kR = min = lim F t [xJ (6.17) 

and every limit of convergent subsequence of x t is a minimum point for F . 

Proof. Let x = (f,g,oo) G Aq such that for all t G (0,1), F t [x\ < 2iiR. 
Then there exists C\ > such that ||V/||l2, ||Vp||i,a||Va;||x,4 < C\. Therefore 
by Poincare's inequality there exists C 2 such that ||/||z,2, ||<7||l2, ||u;|| L 2 < C 2 
and thus by the Banach-Alaoglu theorem the set K = {x G Aq '■ Vt G 
(0, 1), Ft[x\ < 27ri?} is weakly compact. 

1. Let x 4 be a sequence such that min^ F t = F t [x 4 ]. By (I6.12p it follows 
that i*t[xt] < 2tcR and thus by compactness of K there exists a sub- 
sequence x tfc such that x 4fc — *■ x* G Ao. Therefore, by (16. 13)) and the 
above lemma we have that 

minFn < -Fnlx*] < liminf Fu\x. tl 1 = lim min K, = lim min F t . 

Ao fc^o fc fc^O Ao t-w Ao 

2. Fix 5 > and let y G ^4 such that F [y] < min^ F + 5. Then if y,- is 
a sequence satisfying (16. 14)) then 

minFo + 5 > Fo[y] > limsup-Ftfyt] > lim sup min F t = lim min F t . 

Ao t_>o ' i^o A> *-»o Ao 

By the arbitrariness of 5 if follows that 

minFo > lim min F t . 

Ao t-¥0 Ao 



30 



Therefore, by items 1 and 2 it follows that 

min F = lim min F t . 
A t-+0 A 

Furthermore, if x t is sequence as defined in item 1 and x 4fc is a subsequence 
such that x ifc — *■ x* then 

lim min K = minFn < -Fn[x*l < lim min F* 
t-+0 Ao Ao t^O Ao 

which proves that 



min = Fo[x*]. 

Ao 



□ 



Consequently, the small slopes approximation will predict that, for small 
thickness, the configuration of the sheet will converge to a quadratic saddle 
cu = \{u 2 — v 2 ). This is in disagreement with the experimental observation 



that the disks can obtain an arbitrary number of waves [19j. However, if we 



consider the amplitude A n of an n wave configuration, we get the scaling law 

, s 1 / 7T \ 9 TV R 2 , 

A n (R) = -tan — ) R 2 ^ - — , (6.18) 



which does agree with the experiment [19 . 

A key difference between exact isometries and the small slopes approxima- 
tion is captured by the following observation. The curvatures of the surfaces 
given by (16.111) are bounded by max(a,a _1 ) ~ 0(1) independent of R, the 
radius of the disk. This is in contrast with our result (I5.5P which shows that 
for isometric immersions, the curvature has to grow exponentially in the ra- 
dius R. A qualitative expression of this distinction is that the analogs of 
Hilbert's and Efimov's theorems do not hold for the small slopes approxima- 
tion. Consequently, in the next section we look to extend our construction 
by looking for periodic surfaces that are exactly locally isometric to H 2 and 
for small R are well approximated by the small slopes solution. 

7. Periodic Amsler Surfaces 

To mimic the solutions to the small slope approximation we want to con- 
struct a hyperbolic surface that has two asymptotic curves that are straight 
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lines and intersect at the origin of U. Hyperbolic surfaces satisfying this prop 



erty exist, they are called Amsler surfaces [35], and form a one parameter 



family of surfaces Ag that are uniquely determined by the angle 9 between 



the asymptotic lines [25j. As in the construction of surfaces in the small 
slopes approximation, if the angle 9 between the asymptotic lines satisfies 
9 = - we can take the odd periodic extension of the piece of the Amsler sur- 
face bounded between the asymptotic lines and form a periodic profile with 
n waves. We call the surfaces constructed in this manner periodic Amsler 
surfaces. 

The key to generating periodic Amsler surfaces is through the Sine- 
Gordon equation (J3.3P and the similarity transformation z = 2y/uv, (p(z) = 
4>(u, v) which transforms (13. 3B into a Painleve III equation in trigonometric 
form: 

<p"(z) + ^'(z)-sm(cp(z))=0, (7.1) 



where / denotes differentiation with respect to z [25|, |35] . By imposing the 
initial conditions (p(0) = - and <p'(0) = solutions to (17. ip . denoted by <p n (z), 
generate surfaces such that the u and v axis correspond to the asymptotic 
lines of the surface. Therefore, the piece of the surface we are extending 
periodically is parameterized in the first quadrant of the u — v plane and all 
of its geometric quantities are determined by ip n . 

7.1. Maximum radius of periodic Amsler surfaces 

The geodesic equations (13.51) under this similarity transformation take 
the form: 

& + 7k d ^T («*(*>> (f Y - csc(v,> (f) a ) = 



§ + ir, dj t (co%> (f ) 2 - csc(^), (f ) 2 ) = o, (7 ' 2) 



with initial conditions u(0) = v(0) = 0, ^(0) = cos(^), and ^(0) = sin(^) 
for ip G [0, 7r/2]. Notice that these equations have a singularity at the point 
z n = (/?,^ 1 (7r). This is precisely where the immersion becomes singular and 
thus the singular curve is given in the u — v plane by 

v^=^V>. (7.3) 

Once the singular curve is determined we can calculate the maximum ra- 
dius of a periodic Amsler surface with n waves by determining the shortest 
geodesic that intersects the singular curve. 
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' fV> 



Figure 6: Periodic Amsler surfaces with n = 2,3,4 and 5 waves respectively. These 
surfaces are not drawn to scale, but are close to the maximum radii presented in figure 
|9(b)| The coloring of the disks corresponds to contours of arclength data and indicate 
how the different geodesic circles cut from these surfaces would appear. These surfaces 
are in fact not true Amsler surfaces but are discrete Amsler surfaces creating using the 
algorithm presented in [291 ] . 



33 



We can obtain a scaling for ip n 1 (n) by multiplying equation ( 17.1 j) by v 5 '^) 
to get 

1 <%]T = dip d cos if _ 1 2 < c?cos(y?) 

2 dz dz dip z ~ dz 
Integrating, we get that 

* < r m * , , < ,,-h,). (7.5) 



'7r/n a/2(1 - COs(yj)) A/n A/2(cOs(7r/n) - COs(y?) 

Therefore, 

b ( COt ©)-^ (7r) ' (7 ' 6) 

We can use this scaling to prove that we can create arbitrary large 
disks from these periodic Amsler surfaces. First noting that the same ar- 
gument used to obtain (17. 6p implies that lim ip~ 1 (n/2) = oo. Now, let 

n— >oo 

a(t) = x.((u(t),v(t))) be a geodesic lying on a periodic Amsler surface with 
generating angle 9 = n/n. Suppose a(t) starts at the origin and travels to 

the curve y/v = — - — -= — terminating at the point (y,f,Vf) and without loss 

2\/u 

of generality assume that u(t) is a function of v(t). Then, the arclength of 
a(t) satisfies: 




du fduV . f u f 



d(0, («/, v f )) = / ^/l + 2cos(^)— + I — J dv> dv = u f . (7.7) 

Thus, since Uf — > oo as n — > oo we have proved the following proposition. 

Proposition 7.1. Let D be a disk of radius R in the hyperbolic plane. There 
exists a W 2,2 isometric immersion x : D — > U C M. 3 such that U is a subset 
of a periodic Amsler surface. 

Figure [7] is a plot of several solutions to equation (17.11) and includes the 
scaling of the radius of the largest disk with n waves that can be embedded in 
1R 3 . We can see qualitatively that the maximum radius grows logarithmically 
with n, a result that is confirmed by equation (17.61) . The geodesies for the 



particular case when n = 2 are plotted in figure 8(a) along with the geodesic 



circles for n = 3 in figure 8(b) 



34 



Solutions to Painleve Equation 




Figure 7: A plot of ip versus r satisfying the initial conditions (p(0) = tt/u for n = 2, . . . 10. 
The dashed horizontal line corresponds to where (p(r) — n. At this point a principal 
curvature diverges and by choosing a higher value of n a larger disk can be created. The 
inset plot illustrates how the maximum radius scales with n. We can see that the radius 
grows very slowly with n and it looks to be growing approximately at a logarithmic rate. 



Geodesies for an Amsler Surface with Two Waves 
5n 



Geodesic Circles for an Amsler Surface with Three Waves. 





(a) 



(b) 



Figure - 



(a) Plot of the geodesies in the u — v plane of the Amsler surface for n 



The endpoints of the geodesies lie on the curve ypu 



2 waves. 
2 t= , which forms the boundary 
where a principal curvature diverges. 1(b)] Geodesic circles plotted in the u — v plane when 
n = 3. These curves where calculated by contouring the arclength of geodesies. The curve 
in the upper right hand corner is the singular edge of the Amsler surface. 
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7.2. Elastic energy of periodic Amsler surfaces 

To compute the elastic energy using equation (J3.8P we need to determine 
the dependence of u and v on r and \1/. Let j\ and J2 denote unit vectors in 
the u — v plane that are aligned with the u and v axis respectively. If we let 
x denote the parameterization of the Amsler surface, then the pushforward 
or differential of x, denoted x*, is a linear map. Furthermore, since x is an 
isometry and the images of j\ and j 2 under x* make an angle of ir/n we have 
without loss of generality that 

x*ji = ei and x*j 2 = cos(7r/n)ei + sin(7r/n)e 2 , (7.8) 

where e\ and e 2 are the standard basis elements for the tangent plane at 
x(0, 0). Consequently, if we let a(t) = (u(t),v(t)) be a geodesic defined in 
the u — v plane satisfying a'(0) = cos(?/>)Ji + sm (V0J2 then by linearity we 
have that 

x ! „(a / (0)) = (cos(#) + sin(6 l ) cos(7r/n)) e\ + sin(6 l ) sin(7r/n)e2. (7.9) 

Therefore, we can conclude that the polar angle \l/ for this geodesic is given 

by 

/ sin(V') sin(7r/n) \ 

* = arctan — K ^[ l \ l — . (7.10) 

\ cos^^j + sin^j cos{T[/n) J 

Now, we can numerically integrate the elastic bending energy for the 
periodic Amsler surfaces. First, by specifying \l/ and numerically solving 
(17.1 Op for ip we can generate initial conditions for equations (17. 2p which can 
then be numerically solved to determine u(t) and v(t). Then, by fixing r and 
numerically solving the arclength equation 




f du\ , .dudv ( dv\ 

for T we can calculate the values u(T) and v(T) which correspond to the 
coordinates (*&,r). Finally by setting up a mesh on the rectangle (\I/,r) G 
[0, 27r] x [0, R] and using the above process to determine u and v we can 
numerically integrate equation (13. 8p . These results are summarized in figure 

m 
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(a) 



Scaled Bending Energy vs. Radius for Periodic Amsler Surfaces 




(b) 



(a) A representation of the geodesic disks cut out of Amsler surfaces with 
5 waves colored by k\ + k\ in which the radius and polar angle correspond to 



Figure 9: 
n = 2,... 

the geodesic radius r and polar angle \1/. These disks each have a maximum radii of R = 
1.2654, 1.8505, 2.2342, 2.5199 beyond which one of the principal curvatures diverges. We 
can see that these disks have lower energy than there counterparts lying in the hyperboloid 
of revolution. Moreover, the energy is concentrated in small regions near the singular edge 
of the disk. |(b)| The scaled bending energy of periodic Amsler surfaces with n = 2, . . . , 5 
waves plotted versus geodesic radius. The vertical dashed lines correspond to the radius 
where the surface with n waves can no longer be isometrically embedded. 
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8. Discussion 

Free non- Euclidean thin elastic sheets arise in a variety of physical |l|, y, |8[ 
and biological [2|, [ll|, |33[ systems. The morphology of these sheets is usually 
modeled as the equilibrium configurations for an appropriate elastic energy. 
One approach is to model the sheet as an abstract manifold with a prescribed 
target metric g, which is then used to define strains and hence an elastic 
energy for a configuration. 

In this paper we considered the geometric problem of immersing disks 
with constant negative curvature into M 3 . There are two settings for this 
problem which are relevant. One is the problem of "exact" isometric immer- 
sions 

J Dx T Dx = g (8.1) 

where x : D — > M. 3 is the configuration of the center surface. The other 
problem is a small-slopes approximation with an in and out of plane decom- 
position, x = i + eu + e 2 x, g = go + e2 gi with the isometric immersions given 
by 

D X T + D X + Duj T Duj = g 1 . (8.2) 

In the small slopes approximation, we proved that global minimizers of 
the bending energy over solutions of (18. 2p are given by a quadratic saddle 
u(u,v) = \{u 2 — v 2 ). We thus expect that, with decreasing thickness, the 
configuration of the sheet will converge to a two wave solution. This is 
not what is observed experimentally, where the number of waves increases 



with decreasing thickness [19J. Consequently, the physically realized config- 
urations are not given by the global minimizers of the model energy. Note 
however, that the periodic, non-smooth profiles that we construct in Sec. [6] 
do agree qualitatively, and in scaling behavior, with the experimental obser- 
vations. This is a puzzle, since these periodic, non-smooth profiles are not 
global minimizers of the energy, and we need a different selection mechanism 
if these profiles do indeed represent the experimental configurations. 

One way to improve the modeling of the system is to consider the "full" 
Foppl-von Karman energy (12. ip and the associated isometric immersions 
given by (18. ip . For complete surfaces of constant negative curvature (or 
surfaces whose curvature is bounded above by a negative constant), Hilbert's 
theorem [20] (resp. Efimov's theorem 23|) show that there are no analytic 
(resp. C 2 ) solutions of (18. ip . These results are sometimes used as a basis for 
physical intuition that there are no "global" isometric immersions of constant 
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negative curvature surfaces, and this leads to the observed refinement of 



wavelength with decreasing thickness [31 



We investigated this problem and showed (Sec. H] and [5]) there exists 
smooth isometric immersions x : Dr — > R 3 for geodesic disks of arbitrar- 
ily large radius R. Moreover, in elucidating the connection between this 
existence result and the non-existence results above, we provided numeri- 
cal evidence that the maximum principal curvature of such immersions is 
bounded below by a bound which grows exponentially in \/—KqR. 

The existence of smooth, and thus also W 2 ' 2 immersions of arbitrarily 
large disks have consequences for the modeling of free non-Euclidean sheets. 
In particular, the minimizers of the elastic energy (12.1 p must converge as 



t — > to a minimizer of (I2.2p 16l . Il7j , and the energy of the minimizer x* is 

bounded from above by 

£[x*] < t 2 inf B[y] 
ye.4 

where A is the set of all the W 2,2 solutions of (18. ip . This also implies that the 
bending energy of the minimizing configurations should satisfy the inequality 

B[^]<B = inf My] 
ye.4 

for all t > 0, independent of t. This is an apparent contradiction with ex- 
perimental observations where the bending energy of the minimizer scales 
with the thickness t and is well fit by a power law, t _1 , which diverges as 



t — >■ [19[. This suggests that further improvements of the modeling might 
be necessary for a quantitative description of the experimental phenomena. 

The smooth immersions x : Dr — )■ M. 3 (Sees. H] and |5]) do not have the 
rotational ra-fold symmetry (^ — > ^ + 27c/n for the geodesic polar angle 
\l/) of the underlying energy functional, which is also seen in the experi- 
mentally observed configurations. In Sec. [7] We generalize the construction 
of the non-smooth, n-wave small-slopes solutions to obtain isometric im- 
mersions (solutions of (18. ip ) which are also non-smooth and have the same 
symmetry /morphology as the experimentally observed configurations. These 
periodic Amsler surfaces have lower bending energy than the smooth immer- 
sions given by subsets of the pseudosphere and hyperboloids with constant 
negative curvature. 

For each n > 2, there is a radius R n ~ log(n) such that the n-periodic 
Amsler surfaces only exist for a radius < R < R n . This gives a natural 
mechanism for the refinement of the wavelength of the buckling pattern with 
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increasing radius of the disk. However, it does not, at least directly, explain 
the observed refinement with decreasing thickness. 

We have considered the purely geometric model which arises as a limit 
t — y 0. A natural question is the applicability if our results for physical sheets 
where t > 0. Perhaps the refinement of the wavelength with decreasing 
t is an effect which can only be captured by studying the full Foppl-von 
Karman energy (12.11) . For t > 0, we expect that the minimizers of the 
energy will have boundary layers near the edge of the disk and also near the 
singular asymptotic lines, so that the disks satisfy stress and moment balance 
everywhere. The boundary layers near the outer edge of the disk may then 
set the wavelength throughout the disk. A similar phenomenon in stretched 



elastic sheets was studied in this manner in [13 . 

The buckling patterns in non- Euclidean sheets are the result of dynamical 
(albeit slow) processes, and it may be more appropriate to model this type 
of differential growth dynamically, perhaps as a gradient flow for the elastic 
energy. The pattern is then selected for dynamical reasons (fastest growth 
rate, etc.) and not by global minimization of an energy functional. This 
might explain why local but not global extrema for the energy functional 
seem to describe the observed patterns. 

However, neither the boundary layers, nor the gradient flow dynamics, 
can explain how the refinement with decreasing thickness can lead to an 
apparent divergence of the the bending energy £>[x] ~ t~ x . We have shown 
rigorously, for the elastic energy (12. ip . in the limit t — > 0, we should have 
£>[x] — )• Bo < oo. This suggests the need for better models that are carefully 
validated by comparison with experiments in order to obtain a quantitative 
description of non-Euclidean free sheets. 

We conclude our discussion with a short list of mathematical questions 
and avenues for future research that come out of this work - 

1. In Sec. we use an ansatz motivated by numerical experiments to show 
that the solutions to the minimax problem 

0OC = argminmaxcot 2 (0) 

are given by disks cut out of constant curvature hyperboloids of revolu- 
tion, where the admissible set A consists of all C 2 or smoother isometric 
embeddings. We are working on a rigorous (ansatz free) proof of this 
result. 
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2. Our lower bound for the maximum principal curvature (and hence the 
maximum bending energy density) in Sec. can be interpreted as a 
lower bound on the L°° norm of the bending energy density k\ + k\ 
over the admissible set, which consists of isometric immersions. The 
limiting energy and configuration of the non-Euclidean sheet as t — > is 
given by the lower bound for the L l norm of the bending energy density. 
Motivated by this observation, we are investigating lower bounds for 
the L p norms 1 < p < oo to connect the known case p = oo with the 
case of interest p = 1 . 

3. Chebyshev nets give a natural discretization of surfaces with constant 
negative curvature that respects both the intrinsic and the extrinsic ge- 
ometry of the surface. This idea has already been used for constructing 



discrete isometric immersions (i^-surfaces) for such surfaces J29j]. This 
can be interpreted as a t — > limit for the variational problem ( 12. ip . 
We are developing numerical methods for the variational problem given 
by (12. ip with t > that exploits the Chebyshev net structure. 

4. Our results strongly use the fact that the target metric on the sur- 
face has constant negative curvature, which then allows us to naturally 
associate a Chebyshev net with each isometric immersion through an 



asymptotic parameterization [14J. However, we believe that our re- 
sults, in particular the lower bound for the maximum principal curva- 
ture (I5.5p . and the construction of C 1 ' 1 and piecewise smooth isometric 
immersions can also be extended to other target metrics whose curva- 
ture is bounded above by a negative constant. It is natural to study 
these questions, in particular the difference between smooth and C 



i.i 



immersions using the ideas in 44 . 



5. For the periodic Amsler surfaces in Sec. [7J the global structure, i.e. the 
number n of waves is determined by the local structure at the origin. 



The origin is a bifurcation point [45] for the Chebyshev net induced 



by the embedding in the sense that, every point in the complement 
of the origin has precisely two asymptotic directions, while the origin 



has n asymptotic directions with n > 2. The authors of [45[ remark 
that for surfaces whose total positive curvature is bounded by 2n, it is 
always possible to find a global Chebyshev net, except it can have mul- 
tiple bifurcation points. This suggests the following natural questions: 
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Are there non-smooth immersions (C 1 ' 1 or even W 2,2 immersions) with 
multiple bifurcation points? 

The geometry of such surfaces will then serve as a model for the ob- 
served morphology in many non-Euclidean sheets including torn plastic 
[l| and lettuce leaves (lOj which do not have a globally defined number 
of waves, but rather have local buckling behavior which increases the 
number of waves as we approach the boundary. 
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